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Ionization of five inert gas atoms, five alkaHdes, atomic mag- 
nesium, zinc, carbon, nitrogen, and oxygen, and lithium ions by 
electron impact have been considered and the total atomic cross 
section due to single ionization from all shells have been calculated. 
Differential cross sections for ionizaiton of helium and neon by 
electron impact, differential cross section for ionization of helium 
by proton impact, and photoionization of lithium have been calcu- 
lated. Comparison is made with the available experimental data. 
Calculation is based on the Born approximation and a central field 
approximation for the active atomic electron. The active electron 
is described before ejection by a Coulomb function with an effective 
charge given by the Hartree-Fock calculation, and is described by 
a Coulomb function with a charge of unity after ejection. Agreement 
found with a wide range of experimental data for the differential and 
the total cross sections is much better than might be suggested by 
the model of calculation, and in many cases better than the more 
sophisticated calculational models. 

The yield for the production of doubly charged ions in collision 
of electrons with kryton, sodium, potassium, rubidium, cesium, and 
magnesium based on a suggested mechanism has been calculated, 
and in a number of cases satisfactory agreement is found with the 
experimental data. 
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I. INTRODUCTION 


There have been numerous calculations on atomic ionization by charged 
particles. All these calculations are different approximations to a problem 
that cannot be formulated easily. An important case to mention is ionization of 
a hydrogen atom by an electron. Since the quantum mechanical wave function of 
a system of two free electrons and a proton in the asymptotic region of the con- 
figuration space for low impact energies is not known, for the ionization cross 
section no correct functional form or expression which can yield itself to nu- 
merical evaluation is available. At the threshold of ionization the problem has 
been solved classically by Wannier^ with certain validity. 

Here we consider the high impact energy ionization of multiple electron 
atoms by charged particles. Since in this problem the asymptotic form of the 
system's wave function is well known, an expression for the ionization amplitude 
can be found, and with increasing computational labor more and more accurate 
values for the cross section can be obtained. 

A distinction between the ionization of a many electron atom from that of 
the hydrogen atom is the presence of multiple excitation and multiple ionization 
in the many electron case. When resonance scattering for single or multiple 
excitations occurs, giving rise to autoionization, there is an increase in the 
lomc current for certain impact energies, and several maxima are seen in the 
measured ionization cross sections. Another distinction is the interference 
phenomena of the waves that are scattered from different atomic electrons in 
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the case of the many electron atoms. For heavy atoms and molecules whose 
structure resembles a crystal this is an important phenomena. However, dif- 
fraction is more important for elastic scattering. Also, at high impact energies 
where the wave number of the relative motion is many times smaller than the 
atomic radius, diffraction is negligible, and independent particle model as is 
used in this calculation are justifiable. 

Concerning ionization calculation of a many electron atom system, aside 
from the Born approximation, a class of such calculations is the so called 
binary encounter classical calculations. A review of the subject is given by 
Vriens^. In these calculations the target atom is described by a single state. 
Such a model is simpler than the Born approximation where the target atom is 
described by an initial and a final state. A serious defect of a single state 
model IS its inabiHty to produce the important dipole potential field acting on 
the scattered electron which is present in almost all inelastic collisions. This 
potential is due to the fact that quantum mechanically the target is described by 
a mixture of the initial and final states atomic wave functions. Due to this 
potential in accordance with experiment the inelastic cross section at high 
energy behaves as Atii(CE)/E where A and C are some atomic constants and 
E is the impact energy. Since this coupling does not exist in the classical theory, 
this theory is incapable of giving correct results at high impact energies. 
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A correct calculation within the Born approximation requires the knowledge 
of the initial and final atomic state wave functions. While the initial wave func- 
tions are available in the Hartree-Bock approximation, an accurate final state 
wave function, where the ejected electron is not bound to the atom, is not easily 
obtainable. As a result a complete Born calculation does not seem to be 
available. 

McDowell^ and Peach^ have made an approximation to the Born approxima- 
tion by describing the final state as the product of a Hartree-Fock wave function 
for the residual ion and a Coulomb wave function for the ejected electron. 

Implicit m this approximation is that the initial and final states of the atom have 
different Hamiltonians. More specifically, the active atomic electron is in a non 
Coulombic potential before ejection and in a Coulomb potential after ejection. 

They have applied this approximation to the impact ionization calculation of a 
large number of elements. 

e 

A more recent calculation within the Born approximation is due to McGuire 
where the active electron is acted upon by a central Coulomb potential with variable 
charge. The model is used to calculate a number of atomic ionization cross 
sections. 

In the model used here® we describe the active atomic electron before and 
after ejection by a Coulomb wave function. The initial state is then described by 
an hydrogenlike eigenfimction whose quantum numbers are specified by the 
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configuration of the atom and whose Coulomb charge is obtained from a Hartree- 
Fock calculation as will be described later„ This charge from here on will be 
referred to as the effective charge Ze. The final state is described by a Coulomb 
function with central charge unity. Since the effective charge is different from 
unity the two states of the active electron don’t belong to the same Hamiltonian, 
Improvement in the model results if a posteriori the wave functions are made 
orthogonal to each other, as if the two states had the same Hamiltoman. This 
IS done by the conventional method of adding some of the initial state wave func- 
tion to that of the final state. 

An approach similar to ours has been taken by Burhop*^ and a number of 
other authors. According to these authors the active electron is both initially 
and finally in a Coulomb field with the same effective charge; thus no orthogo- 
nalization of the states is necessary. We believe that this is not a good approx- 
imation in view of the fact that the long range dipole potential mentioned 
previously plays an important role at high impact energies, and in the model 
of Burhop and others the dipole moment in the potential is smaller than it 
should be. To see this we notice that the time spent by the atomic electron in 
traversing the diameter of the atom before ejection is much smaller than the 
time spent by the scattered electron in the dipole potential field. An appropriate 
final state for the atomic electron would then be a Coulomb function with central 
charge equal to unity. The effective charges used by Burhop and others are 
greater than unity. This would lead to a dipole moment for the long range potential 
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smaller than the dipole moment used m the present paper, and this leads to cross 
sectional values which at high energies will be smaller than those calculated 
here. Comparison with the experimental data to be discussed later shows that we 
have satisfactory agreements with measurements for a large number of cases at 
high impact energies. 

We now summarize the results obtained here. Among the inert gases the 
ionization cross section by electron for helium, neon, argon, krypton, and xenon 
have been calculated. In alkalides similar calculation has been done for lithium, 
sodium, potassium, rubidium, and cesium. Total cross section for 
electron Impact ionization of magnesium, zinc, carbon, oxygen, and nitrogen 
have also been calculated and given. In all cases the ionization cross section of 
not only the outermost shells but all the shells that have contributed appreciably 
to the cross section have been calculated. All calculated results have been 
compared with those of experiments and interpretations have been made accordingly. 

A critical test for the validity of the model is comparisons between calcu- 
lated and measured differential cross sections per unit energy range of the 
ejected electrons. Here this comparison is made for protons on helium, and 
electrons on helium and neon where the measurements have recently become 
available. 

Similarly, it is instructive to compare calculated photoionization cross 
sections based on the present model with measurement. Since photoionization 
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cross section has the same functional form as the high impact energy cross 
section per unit energy range of the ejected electrons, this comparison will 
also give some information about the accuracy of the model. Here we have cal- 
culated photoionization of lithium, and have made a comparison with the experi- 
mental results. 

It will also be interesting to see the degree of validity of the model for ioni- 
zation of ions. For this reason we have calculated ionization of lithium ion by 
electron impact and have compared calculation with measurement. 

Finally, we like to mention the limitations of the present model, and also of 
the Born approximation. Since the difference between the two Hamiltonians of 
the atomic states in the present model is proportional to the difference in the 
two Coulombic charges, when the initial effective charge is greatly different 
from unity, the approximation will not be a good one. Then the model provides 
a better approximation for the outer shell electron lomzation than the inner 
shells. From the point of view of the relativistic theory, for some inner shells 
Z e/137 is not small compared to unity, and the Dirac's theory should be used to 
obtain the atomic wave functions. On the other hand the central field approxima- 
tion is invalid for description of the outer shell electron wave functions in the 
heavy elements. Some other approximations, like that of Thomas Fermi model, 
should be used to obtain these wave functions. 
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n. FORMUIATION 


The cross section per unit energy range of the ejected electrons, do-/de, in 
the Born approximation and according to the model discussed before in ionizing 
collision of a charged particle with charge Z'e, e being the electronic charge, 
with an atom with an active electron specified by the principal and the azimuthal 
quantum numbers n and ^ is given by 



with cr the total cross section in umts of a^, e the energy of the ejected electron, 
E and M the relative energy and reduced mass of the colliding system, and 
the electronic mass, kj is the wave number of the relative motion before col- 
lision related to E by 

and kj is the wave number of the relative motion after collision related to k^ by 

/?,*■- (S) 

with AE the excitation energy. Tiq (k^ “ ^2) momentum transfer be- 

tween the particles, m is the absolute value of the magnetic quantum number of 
the atomic electron, and k is a unit vector in the direction of the ejected electron. 
The final state vector | f ^ is given by 

1 "iy. j K )> t y i'y) 
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where ]k^ is a Coulomb wave function with charge imity normalized such that 
asymptotically it behaves as ejcp (ik • r) with r the position vector of 

the atomic electron. The constant y is given by 




(5) 


so that 0, as required. The momentum of the ejected electron is k 

and is related to e by the relation 


/k"/ = 


( 6 ) 


It is easy to see thato- has the correct dimension as required. 


From (4) and (5) we have that 

</ / - </T| 

— {^1-^ {-'’/'»> je.'" ^ I ( 7 ) 

It is easier to evaluate the first matrix on the right hand side of (7) in parabolic 
coordinates and then transform it to the spherical coordinates. If the state 
vectors m the parabolic coordinates be specified by Innjtn^, with the usual 
quantum number in parabolic coordinates, then it follows immediately that 



where the coefficients ^nn^m are given through the Wigner's 3j symbol by 
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witlL n 2 other quantum number in parabolic coordinates related to by the 
relation + n2 = n- m-l. The value of the matrix j exp ( i q ’ r ) | nn^ is 
given below. 
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In these equations a. = Z^ Z^ being the effective charge, j6 = Z/ka^, where 
a^ is the Bohr radius and Z = 1 is the charge of the Coulomb field of the ejected 
electron, and A, B, and D are defined by 




(13) 


The constants C (,u. 4 » M 3 ) are given elsewhere. 4>^ is the angle of azimuth of 
q with respect to k as the z-axis. The lower limits of all the demi variable 
integers in ( 11 ) are zero. 


Bor evaluation of the second term on the right hand side of (7) we notice that 
I n-tm^can be obtained by malcmg | q | in the first terra equal to zero. The 
second factor in this term can be evaluated by the usual methods; although the 
calculation is lengthy in algebra for large n values used here. The matrices in 
(7) involving k have been evaluated by taking the z-axis along k. Then the second 
term on the right hand side of (7) vanishes unless m == 0. 

In this way the matrix in (1) is evaluated and dcr/d( e/Ry) found accordingly. 
An integration of do-/d(e/Ry) with respect to e/Ry will yield the total ionization 
cross section for ejection of an electron in a given atomic shell. 
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The effective charge for the atomic electron before ejection is fonnd in the 

following way. Let Z be the atomic number of the atom, then the effective charge 

is given by = Z - S, where S is the screening parameter representing 

screening of the nucleus by other atomic electrons. S is defined by Hartree^^ 

through the relation R = ^ some linear scale of the atomic wave 

function under consideration, and the corresponding value for the atomic 

hydrogen. It will be a good approximation to take R to be equal to the mean 

12 

atomic radius r . Froese has assumed this to be the case, and gives S for 

several elements. We have used here values of S given by this author. For 

13 

ionization of ions we have used the similar parameters given by Naqvi. 

Effective charges and the ionization potentials of different atomic shells for the 
elements used here are given in Appendix I. 

I 

An examination of Eqs. (10) through (13) shows that the integrand on the 
right hand side of (1) depends only on q, k, and the angle between q and k . 

Then for evaluation of (1) a two dimentional numerical integration, namely inte- 
gration with respect to this angle and q, is necessary. An additional integration 
with respect to e will give the total cross section. 

In some measurements cross section per unit energy range of the ejected 
electron and per unit solid angle of the scattered electron is desirable. We give 
an expression for this quantity for the future use. Remembering the relationship 
between q and we have that qdq = k.^ 2 ’^ 0 .^/{ 2 Tr), where dRg is the element of 
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the solid angle of the scattered particle with respect to the incident particle as 
the z-axis. From (1) the desired cross section is given by 


r-/4. 



where the two dimensional angular integration is with respect to as the 
z-axis. While in (1) the integral is cyclic with respect to the angle of azimuth 
of k , this is not the case here. If and represent angles of azimuth of 
and k with respect to k^^ as the z-axis, the matrix in (14) contains terms with 
argument (0^ " This indicates that the left hand side of (14) is independent 
of 


For ionization of ions we have used plane wave for the incident, and Coulomb 
wave for the ejected, electrons. The use of the plane wave instead of the Coulomb 
wave for the incident particle is justifiable at high impact energies and the error 
incurred is of the same order as the error due to the use of the Born approxi- 
mation, But the ejected electron wave function should always be described by a 
Coulomb wave. 


For evaluation of the photoionization cross section we have used the well 
known formula^"^ 
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in winch v is the frequency of the incident radiation, k the wave number of the 
ejected electron related to v by the conservation of energy through k^ + 

I = hv, I being the ionization potential of the ejected atomic electron, | C} is the 
eigenstate of the atomic electron and its degeneracy number. [ is the 
Coulomb function of the atomic electron with normalization and effective charge 
described previously, k is a unit vector m the direction of k. It is seen that o-^ 
has the dimension of length squared. 

A relationship exists between photo ionization cross section and the dipole 
approximation of the charged particle impact ionization. This can be obtained 
by expandmg exp (i q • r) in (1), retaining the first non-vanishing term of the 
matrix, and employing (15) . Then the differential cross section in the dipole 
approximation will be related to <t„ through 

dcr^ _ Z 

E(r-,e) 

where a ~ 1/137 is the fine structure constant of the Coulomb potential, and I is 
the ionization potential. 

For small momentum transfer where the dipole approximation is valid we 
can write + k 2 = 2kj, and - kj = AE/(2a^ k^Ry), where AE = I + e is the 
excitation energy. In this way (16) reduces t 

— ^ 

' < 17 ) 




X, 


-f 


(16) 
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giving the differential cross section in terms of the photoxonization cross section 
o-^. Obviously v = A E/h. In employing (17) note should be taken of the fact that 
cr^ should correspond to low ejected electron energies, i.e., energies of the order 
of several threshold units, since the dipole approximation is valid for distant 
collisions, or low energy ejected electrons. 

in. RESULTS AND DISCUSSION 

The results of our calculations are presented in the following 23 figures. 
They are compared with the experimental data and with the results of some 
other representative quantum mechanical calculations. 

The experimental data up to 1966 have been collected and reviewed in an 

15 

article by Kieffer and Dunn. Much of the experimental data used here has 
been taken from this review and referred to this reference accordingly. 


In comparing the results of calculations with the experimental data note 
should be taken of the fact that in almost all measurements the quantity meas- 
ured is the amount of ions produced. The measured ionic current I^ is given by 




-f- 


(18) 


where is the current due to the atoms which have lost n electrons. Some 
of the current is produced by direct collision and some is due to after collision 
processes such as autoionization or Augei' effect. 


15 



There has been a number of measurements in which ions with specific charges 
on them, like singly or doubly charged ions, have been collected and measured. 

For all experiments mentioned in this paper all these measurements are relative 
while the measurement of the total ionic current given by (18) is absolute. Also, 
these measurements have less accuracy than the measurement of the total ionic 
current. 

It then becomes obvious that in comparing calculation with measurement 
multiple ionization as well as autoionization should be taken into account. In 
the calculation reported here we have limited ourselves to single ionization 
from the individual shells only. Since the single ionization is the dominant 
process throughout the impact energy range, satisfactory agreement with ex- 
perimental data is obtained in many cases. As an example, from the experi- 

r 

mental evidence the production of singly charged ions from all shells at high 
impact energies accoimts for 98% in helium, almost all in lithium, and 90% 
in neon, of the total ionic current. Then, for light elements up to neon we 
justifiably can compare our calculated cross section with measured total ioni- 
zation cross section. For heavier elements this comparison is less justifiable 
due to the presence of many autoionization states. 

In production of doubly charged ions four processes can take place. Three 
of these are double ionizations, simultaneous ionization and excitation, and double 
excitations. The latter two in certain cases may lead to the ejection of one for 
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the first case, and two for the second case, of additional electrons. The fourth 
process is a single inner shell ionization with subsequent cascading of the outer 
shell electrons and ejection of one or more electrons, called Auger electrons. 
Due to the large cross section of the single ionization compared to the multiple 
ionization, the fourth process at high impact energy is usually the dominant 
process compared to the other three. This is supported by the fact that our cal- 
culated cross section for single ionization of some inner shells almost equals the 
measured cross section for production of doubly charged ions in many cases. 

We now give a description of the results. 

A. Differential Cross Sections 

Before giving a description of the total cross sections it is instructive and 
a critical test for the theory to present results for the calculated differential 
cross sections and a comparison with the experimental data available. By the 
differential cross section we mean here cross section per unit energy range of 
the ejected electrons. It will be a better test for the theory to compare cross 
sections per unit energy range of the ejected electron and per unit solid angle 
of the scattered electron with the experimental data, but this will not be done 
here. For the future use the formula for this cross section is given by (14). 

In Fig. 1 the calculated differential cross sections for ionization of helium 
and neon by electron impact are given and compared with the recent measure- 
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ments of Opal, Peterson, and Beaty. These are relative measurements 
normalized by the experimenters . 

In the case of helium our calculation agrees fairly well with the experi- 
mental data except for the up turn seen at 50 eV m the experimental curve 
corresponding to the 100 eV primary energy. This up turn is most likelj^ due 
to the exchange effect not taken into account m our calculation and which is 
important at the low impact energies. 

In the ease of neon the agreement is not as satisfactory as in the case of 
helium. In particular, the calculated curve seems to be flat in the region of low 
ejected electron energy. This is characteristic of the p-shell electron ionization 
in our model whose cross section dominates the neon, and also the carbon, nitro- 
gen, and oxygen lomzation differential cross sections to be discussed shortly. 

In Fig. 2 the calculated differential cross section for ionization of helium 

by proton impact is given and compared with the measurements of Rudd, Sautter, 
17 

and Bailey. As can be seen the agreement is quite satisfactory. A more elab- 
orate calculation on differential cross section for the ionization of helium by 

electron and proton impact with an accurate wave function has been done by Bell, 

18 

Freeston, and Kingston. The amount of agreement of the results of Bell et al 
with the data of Rudd et al is comparable to the agreement that we have found in 
our comparison with the data of Rudd et al. This indicates that a very accurate 
wave function is not necessary to calculate this differential cross section. 
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A similar test for the accuracy of the model employed here is a comparison 
of the measured and calculated photoionization cross sections. Using our data 

3 

this comparison has been made for lithium by McDoweU. Here this comparison 

is reproduced and is shown in Fig. 3. The experimental data belongs to Hudson 
19 

and Carter. The calculated cross section is lower for low, and is higher for 
high, ejected electron energy compared to the experimental data. 


B. Inert Gases 

Our results for the inert gases are given in Figs. 4-9. For these gases we 

restrict ourselves chiefly to the calculation for production of singly charged ions. 

20 

Source for ail the experimental data, except those of Goudin and Hagemann, and 
21 

Schram, is reference 15. It is seen that about the peak of the ionization curve, 
the experimental data are not consistent with each other, and they have consid- 
erable spread. 


In the cases of helium, neon, and argon the experimental data displayed 
represent total experimental cross sections and are respectively about 2%, 10% 
and 12-15% higher than the observed single ionization cross sections in the 
energy range from 500 to several thousand eV. 

In Figures 4, 5, 6 we have presented the results of our calculation and cal- 
4 22 

culation. of Peach. ’ For neon and argon our calculation clearly gives better 


agreement with experimental data than those of Peach, 
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For krypton and xenon (Figs. 7, 8, and 9) the data of Tate and Smith are 

normalized to the multiple ion production measurements of Rapp and Golden at 

21 

400-500 and 500 eV, respectively (Cf. Ref. 15). Similarly, the data of Schram 
are normalized to an earlier work of Schram, de Heer, Wiel and Kistemaker 
(Cf. Ref. 15). 

For krypton it is seen that in their energy range the data of Schram comple- 
ments those of Tate et al. T(4s + 4p) is our calculated curve, and as is seen at high 
impact energies, it is in reasonable agreement with the data of Schram. 

For xenon (Pig. 8) the two measurements are in some disagreement with 
each other. Our calculation agrees better with measurement of Schram, and is 
somewhat higher than these measurements. 

In Fig. 9 we compare our results for a single 3d shell ionization in krypton 

23 21 

with the data of Tate and Smith and Schram for production of doubly charged 
ions in this element. We argue that the mam source of the doubly charged ions 
at high impact energy is a single 3d shell ionization with a subsequent ejection 
of an Auger electron. Our results agree quite well with the high energy data of 
Schram. The disagreement with the data of Tate and Smith may be due to the 
production of doubly charged ions due to other processes important at low energies. 

C. Alkalides 

In Figures 10-14 ionization cross sections of allcalides ranging from lithium 
to cesium have been given. In Fig. 10 our results for both shell ionization clearly 
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gives better agreement with the absolute measurement of McFarland and Kinny 

28 

than the results of McDowell at al marked MMP(2s). At high impact energies 
both calculated results seem to be lower than the results of measurements. 

In Fig. 11 cross sections for production of singly and doubly charged sodium 
ions are presented. In accordance to what was said before production of the Na^^ 
ions is mainly due to the single ionization of the 2s shell and a subsequent Auger 
process. This is confirmed by the fact that a- (2s) which gives cross section for 
the ejection of a 2s electron shown in the figure by T(2s) is in agreement with 
cross section for the production of doubly charged ions. 

In Figs. 12-14 ionization of potassium, rubidium, and cesium are presented. 
T(3s), T(4s), and T(5s) in the three figures represent twice the calculated cross 
sections for single ionization of the inner s-shell electrons, since these single 
ionizations give rise to doubly charged ions. In these figures the experimental 
curves for the production of doubly charged ions have lower thresholds compared 
to T(3s), T(4s), and T(5s) curves. Also, they give larger cross sections at low 
impact energies. This can partially be explained by the importance of autoioni- 
zation at low energies. In particular, autoionizations permit formation of doubly 
charged ions below the threshold of the inner s shell electrons. At higher energies 
autoionization is probably not as important and the Coster-Kronig type Auger 
electrons can account for the majority of the doubly charged ions. 
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It should be emphasized that the experimental curves m Pigs. 10-14, unless 


otherwise specified, are for production of the total ionic current. While the pro- 
duction of ions with charge greater than two is negligible for Li, Na, and K, Tate 
and Smith have shown that at 500 eV 13% for rubidium and 18% for cesium of 
the total ionic current are due to triply or higher charged ions. Thus in compar- 
ing the experimental data of McParland and Kinney with our total cross sections 
T(4s +4p + 5s) for Rb and T(5s + 5p + 6s) for Cs an appropriate reduction in their 
cross sections should be made, but this has not been done here. 

It will be interesting to notice how the results of McFarland and Kinney, 
which are most recent, absolute, and are extended to high impact energies, 
compare to our results for the five alkalides given in Figs. 10-14. It is noticed 
that for Li, Na, and Cs the experimental values lie above, and for K and Rb they 
lie below, the theoretical results. Although we have employed a crude model for 
calculation, no consistent discrepancy is found with measurements. More refined 
calculation should give better agreement with measurements. 

There seems difficulty in reconciling our results with the data of Tate and 
Smith for high energy ionization of Rb and Cs (Cf. Figs. 13 and 14) as is seen in 
the "flattening" of the curves due to these data. Similarly, as is seen in Figs. 13 
and 14, while the T(5s) curve representing single 5s ionization correctly gives 
the production of doubly charged cesium ions, the T(4s) curve does not give the 
production of the doubly charged rubidium ions. To resolve the differences, and 
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since the ineasurement of Tate aad Stnith, is 36 yeaxs old, a new measurement 
for the relevant cases is desirable. 

The irregular variation in the ion production in K, Rb, and Cs at low impact 
energy is due to numerous autoionizing states in these atoms which has not been 
taken into account in our calculations. 

D. Magnesium and Zinc 

In Figures 15 and 16 our results for ionization of magnesium and zinc are 
shown. As is seen in Fig. 15 our calculated curve marked T(3s) for single ion- 
ization of magnesium is in fair agreement with the measurement at high impact 

energies, but near the peak calculated results are substantially lower than the 

34 

measurement. Measurements of Kaneko {not shown in the figure) indicate the 
presence of a number of autoionization lines from threshold up to the peak of the 
ionization curve. 

In Fig. 15 it is also seen that the cross section for production of doubly 
charged ions compared to singly charged ions in magnesium is much larger than 
the same quantity for sodium. While for sodium the doubly charged current 
accounts for only 10% of the total cross section at high emergies, in magnesium 
it is 50%. The difference is most likely due to the fact that in Mg the ejection of 
a 2p electron can be followed by the emission of a 3s electron through the Auger 
process. This cannot take place m sodium since the energy gap between 2p and 
3s in. sodium is not wide enough to make the process operational, and also due 
to the fact that there is only one 3s electron in sodium. 
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In Fig. 15 the curve Mg 2 + lies higher than the calculated curve T(2s +2p) 

which is assumed to give the cross section for the Auger process. The discrep- 
ancy may be due to a number of reasons and cannot be explained easily. 

Figure 16 gives calculated cross sections for ejection of 4s, 3d, and 3p 

electrons in zinc. No measurement on ionization of zinc is available. Photo- 

35 

ionization data on Zinc indicate that the element is rich in auto ionization states 
at the low energy side of the peak of the ionization cross section. We then 
anticipate that the peak in the measured ionic current will be higher than the 
peak of the curve marked T(3p + 3d+4s). As usual we have assumed that 3p shell 
electron ionization will lead to the production of doubly charged ions through the 
Auger process. For this reason <^(3p) has been multiplied by a factor of 2 before 
being added to the cross section for the total ionic current marked by T(3p + 3d +4s). 

36 

The shape of the curve found by Fiquet-Fayard et al for production of Zn’^'®’ 
(not shown m the figure) is similar to T{3p) above 500 eV. Since the cross sec- 
tion as a function of energy at high impact energy according to the T(3p) curve is 
given as E“^(AlnE +B) with A and B constants, the similarity suggests that Zn^ 
production at high energy is dominated by the Auger process due to single ioni- 
zation. At lower energies other processes should be taken into account. This 
would probably make the curve for the production of Zn to lie higher than T(3p) 
at low energies. 
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Because of the small size of the Biquet-Fayard figures, their reproduction 
is difficult and is not done here, 

E. Carbon, Nitrogen, and Oxygen 

Figs. 17-19 gives our results for single ionization of carbon, nitrogen and 
oxygen. No experimental data is available for carbon. For nitrogen and oxygen 
the experimental data are in as much disagreement vdth each other as are dif- 
ferent calculations. 

It is important to note that in these three elements contribution to the ioniza- 
tion cross section of the 2s shell is almost as important as the 2p shell, and in 
any calculation th^ two shells should be taken into account. The earliest calcu- 
lation due to Seaton, S(2p) in the figure, contains contribution of the 2p shell only. 

In Figs. 20-22 cross sections per unit energy range of the ejected electrons 
for electron impact ionization of carbon, nitrogen, and oxygen are given. Because 
of the abundance of these elements in the atmosphere and other parts of the cos- 
mos, these cross sections will be useful for a number of problems related to the 
stopping power of the atmospheric and cosmic gases. The agreement found in 
Fig. 1 indicates the likelihood of the accuracy of the differential cross sections 
given m Figures 20-22. 


F. Ionization of Ions 

In the ionization of ions the incident particle should be described by a Coulomb 
wave function. At high impact energies the Coulomb wave can be replaced by a 
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plaae wave. It can be shown that the resultant inaccuracy is comparable to the 

inaccuracy in the Born approximation. This has been verified by Peart, Walton, 

38 

and Bolder who have compared their absolute cross section measurement of 
electron impact ionization of He'*’ with a Born calculation where a plane wave is 
used for the incident electron. 

Here we apply our model to the ionization of a non-hydrogenlike ion, namely 
Li'*'. Effective charge and lomzation potential of Li'*' necessary in our calculation 
is given in the Appendix. 

In Figure' 23 we have constructed a Bethe plot for ionization of Li by electron 

impact. Measured cross sections are by Peart et al (opt. cit.) , and by Peart and 
39 

Bolder. The measurement is well extended into the relativistic region. When 

the projectile has relativistic energies, the target, excluding the inner shell of 

heavy atoms and ions, can still be treated non-reiativistically. The relativistic 

40 

correction for the projectile has been made by Bethe. The non-re lativistie 

41 

high energy inelastic atomic cross section according to Bethe is given by 
cr = (A in E + B) with E the relative kinetic energy of the system, and A and 
B some atomic constants. For relativistic projectiles this expression should be 
replaced by 

( 19 ) 
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with m and v the rest mass and velocity of the projectile, ji = v/c, c being the 
speed of light, and A and B the non-relativistic atomic constants* 

In Fig. 23 the straight line R is a plot of the left hand side versus the ex- 
pression in the curly bracket in (19) . The slope and the intersection of R gives 
A and B. It has been obtained by constructing a segment of R in the non- 
relativistic high energy region, using for cr the computed Born cross section. 

By extension of the segment R has been obtained. The curve NR has been ob- 
tained by using for o- the computed Born cross section throughout the range of 
the curve. The departure of NR from R curve at high impact energy is a measure 
of the amount of the relativistic correction. 

If the relativistic correction were not necessary the experimental points 
should fall on NR curve, assuming the accuracy of our non-relativistic calcula- 
tion. With relativistic correction the esqperimental data should agree with R. 
Although the experimental data are extended to the relativistic region because 
of their error bars the correctness of the relativistic correction cannot be 
verified, and this verification is left to the future. 

Before closing this section we make three remarks. It is of interest to know 
the relative contribution of different atomic shells to the ionization cross section. 
There is no simple quantum mechanical expression for this relative contribution, 
but an order of magnitude relation can be obtained using the Thomson's classical 
formula for ionization. According to this formula the ratio of the cross sections 
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Q' and Q ^ for two elements 1 and 2 evaluated at equal tlir eshold units of ener- 
gies IS given by ^ 

where and I 2 are the ionization potentials of the two elements. Thus, as the 
ionization potentials of the inner shells get larger, lomzation cross section of 
these shells get smaller, approximately as the inverse square of the ioniza- 
tion potentials. Contribution of the inner shell ionization in many cases can then 
be neglected. 

A second remark is about construction of the Bethe plot for many electron 
atom ionization. Since high energy non-relativistic cross section is given by 
o- = E (A In E -r B), a plot, called the Bethe plot, of Ecr versus In E falls on a 
straight line. Eor heavy elements when the incident energy has reached many 
times the ionization potential of the outer shell electrons, corresponding to the 
validity of the above expression, the energy may not have reached' the threshold 
energy of the inner shell electrons. Thus construction of Bethe curves for heavy 
elements in the energy range of interest is not possible. For lighter elements 
this is possible, and the two parameters of the curve give all the necessary 
information for the high impact energy total cross sections. 

Finally, we have established that Auger transition in production of doubly 
charged ions is a dominating process. This is supported by the present calcu- 
lation which shows that for several elements the inner shell ionization responsible 
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for the Auger transition yields cross sections very close to that observed for the 

21 44 

production of doubly charged ions. In addition, Schrano et al ’ and Piquet- 
36 

Payard et al have shown that cross section for this production behaves as 
E~^(A InE +B) with A and B comparable to each other in their magnitudes. Since 
in multiple ionization or excitation A is very small compared to B, our assump- 
tion is verified. The cases of K and Rb may however be exceptions as the meas- 
urement of Fiquet-Fayard et al above 500 eV indicates that the cross section for 
production of and Rb"’"'*' (not shown m the figures) falls off faster as a function 
of energy than the curves T(3s) and T(4s) in Pigs. 12 and 13 indicate. The dis- 
crepancy may be due to a number of reasons which we will not discuss here. 

IV. CONCLUSION 

1. By the use of a simple model we have calculated a number of ionization 
cross sections in fair agreement with measurements. Its simplicity and relative 
accuracy makes it a useful calculational tool for further ionization calculations. 
We also have shown that the model can be used for ionization of ions with an 
accuracy comparable to the atomic ionization. 

2. In comparing our results with the experimental data for differential 
ionization cross sections we have found that agreement is comparable to the 
agreement of the more sophisticated calculations with accurate wave functions 
We conclude that the accuracy of the differential ionization cross section is not 
very sensitive to the accuracy of the wave functions. 
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3. We have shown that in a number of elements production of doubly charged 
ions is mainly due to Auger processes whose cross section is given by certain 
inner shell ionization cross sections. 
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FIGURE CAPTIONS 


Fig. 1. 

Fig. 2. 

Fig. 3. 
Fig. 4. 

Fig. 5. 
Fig. 6. 


Cross section per unit energy range of ejected electrons in ionization 
of helium and neon by electron impact. The energy of the primary in 
eV is marked on each curve. The measurements are due to Opal, 
Peterson, and Beaty (Ref. 16). 

Cross section per unit energy range of ejected electrons in ionization 
of helium by proton impact. The experimental data are due to the 
absolute measurement of Rudd, Sautter, and Bailey (Ref. 17), 

Photoionization of lithium. Measurement of Hudson and Carter 
(Ref. 19) is compared with our calculation. 

Ionization of helium by electron impact. In this and all the following 
figures through Fig. 19, and Fig. 22, a is the ionization cross section 
in units of and E is the impact energy in units of eV. All the 
experimental data in this figure and the following two figures are 
absolute measurements. For a description see Ref. 15, except the 
data of Gaudin and Hagemann which is described in Ref. 20 . The two 
solid lines are calculations of Peach (Ref. 22), and ours. 

Ionization of neon by electron impact. See Refs. 15 and 20 for a de- 
scription of the measurements. T(2p), T(2s), and T(2s + 2p) represent 
present calculations and give cr(2p), cr(2s), and the sum of o-(2p) and 
cr(2s), respectively. P(2s + 2p) is calculation of Peach (Ref, 4) and 
gives the sum of cr(2p) ando-(2s) plus autoionization. 

Ionization of argon by electron impact. See Refs. 15 and 20 for a 
description of the measurements. T(3p), T(3s), and T(3s + 3p) repre- 
sent the present calculations and they give cr(3p), a-(3s) and the sum 
ofcr(3p) and(r(3s) respectively. P(3p) andP(3s) are calculations of 
Peach (Ref, 4) and they give cr(3p) and cr(3s) , respectively. In P(3s) 
autoionization has also been taken into account. 
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Pig. 7. Single ionization of krypton by electron impact. Kr curve gives the 

measured cross section for production of Kr'*'. The relative measure- 
ments of Tate and Smith (Ref. 23), and Schram (Ref. 21) are normal- 
ized as described in the text. T(4p), T(4s), and T(4s + 4p) represent 
the present calculations and they giveo-(4p), cr(4s), and the sum of 
o-(4p) and cr(4s), respectively. 

Fig. 8. Single ionization of xenon by electron impact. Xe'^ gives relative 

measured cross section for production of singly charged xenon ions. 

See Fig. 7 for references. T{5p), T(5s), andT(5s + 5p) represent the 
present calculations and give cr(5p) , o^(5s) and the sum of o(5p) and 
cr(5s) , respectively. 

Fig. 9. Production of doubly charged krypton ions by electron impact. 

curve gives twice the relative measured cross section for 'production 
of doubly charged ions, see Fig. 7 for references. Present calculation 

t 

gives 2a-(3d) with cr(3d) the cross section for single ionization of a 3d 
shell electron. The colliding atom becomes doubly ionized by sub- 
sequent ejection for an Auger electron. 

Fig. 10. Ionization of lithium by electron impact. Measurements of McFarland 
and Kinney (Refs. 15, 24, 25), and Zapesochnyi and AleksakMn (Ref. 26) 
are absolute total eiqperimental cross sections; those of Brink (Ref. 27) 
are relative normalized to the results of McFarland et al. The two 
curves marked T are the present calculations; one for ionization of 
the outer shell, and the other for ionization of both inner and outer 
shells of lithium. MMP(2s) is the outer shell calculation of McDowell 
et al (Refs. 3, 28). 

Fig. 11. Ionization of sodium by electron impact. Measurements of McFarland 
and Kinney (Refs. 15, 24, 25), and Zapesochnyi and Aleksakhin (Ref. 26) 
are absolute total experimental cross sections. Those of Tate and 
Smith (Ref. 23), and Brink (Ref. 27) are the cross sections for the 
production of the ions Na"^ and They are relative and are normal- 

ized at 500 eV to the results of McFarland and Kinney. Curves marked 
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by T represent the present calculation. T(3s) gives ct-(3s), 3s being 
the outermost shell, T(2s) gives 2cr(2s), andT(2p + 3s) gives the sum 
of cr(3s) and o-(2p) and is to be compared to the experimental Na"^ 
curves. T(2s + 2p+3s) gives the sum of cr(3s), cr(2p) and2v(2s), and 
must be compared with the measured total ionization cross section. 
For an explanation see the text. BBP(3s) is calculation of Bates et al 
(Ref. 29) for cr{3s), andP(2p + 3s) is calculation of Peach (Ref. 4) for 
the sum of o-(3s) and cr(2p). 

Fig. 12. Ionization of potassium by electron impact. The measurement of 

Korchevoi and Prozonski (Ref. 30) is absolute. For other measure- 
ments see Fig. 11. T curves represent the present calculations, 

T(4s), T(3p) and T(3s) gives cr(4s) and o-(3p), and 2cr(3s), respectively. 
T(3p + 4s) gives the sum of cr(4s) and a(3p), and is to be compared 
with the experimental curves for K*", and T(3s + 3p+4s) gives the sum 
of cr(4s) , cr(3p), and 2o-(3s) and is to be compared with the total ex- 
perimental cross sections of McFarland and Kinney. 

Fig. 13. Ionization of rubidium by electron impact. For explanation of meas- 
urements see Fig. 12, T curves represent the present calculations. 
T(5s), T(4p), and T(4s) give ct(5s), cr(4p), and 2cr(4s), respectively. 
T(4p-i-5s) gives the sum of o-(5s) and o-(4p), and is to be compared 
with the Rb"^ curve of Tate and Smith. T(4s + 4p+ 5s) gives the sum 
of o-(5s) , cr(4p), and 2cr(4s) but in this case it should be smaller than 
the total experimental cross sections of McFarland and Kinney, see 
text. 

Fig. 14. Ionization of cesium by electron impact. For measurement of Heil 

and Scott see Ref. 31. For other measurements see Fig. 12. T curves 
represent the present calculations. T(6s), T(5p), and T(5s) give cr(6s), 
cr(5p) , and 2o-(5s) , respectively. T(5p h- 6s) gives the sum of cr(6s) and 
<r(5p) and should be compared with the measurement for Cs^. 

T(5s + 5p+6s) gives the sum of cr(6s), o-(5p), and 2o-(5s) and should be 
compared to the total experimental cross sections of McFarland and 
Kinney, see text. 
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Fig. 15. Ionization of magnesium by electron impact. Mg"^ and curves 

gives cross sections for production of singly and doubly charged ions; 
cross section for the latter curve being multiplied by a factor of 20. 

The measurement of Okuno et al (Ref. 32) is the absolute total ioniza- 
tion cross section. The relative measurements of Okudair-a et al 
{Ref. 33) were normalized by the experimenters at 500 eV to the re- 
sults of Okuno et al. T curves represent the present calculations. 

T(3s), T(2s + 2p), and T(2s + 2p+ 3s) give a(3s), 20 [cr(2p) + cr(2s)] , and 
o^3s) + 2cr(2p) + 2cr(2s), respectively. P(3s) andP(2p) are calculations 
of Peach (Ref. 4). They give cr(3s) and 20cr(2p), respectively. In 
P(2p) curve the auto ionization has also been taken into account, 
accounting for its dip. 

Fig. 16. Ionization of zinc by electron impact. T curves represent the present 
calculations. T(4s), T(3d), and T(3p) give cr(4s), cr(3d), and40o-(3p), 
respectively. T(3d + 4s) is the sum of cr(4s) and a(3d), and T(3p+3d + 4s) 
is the, sum of cr(4s) , o-(3d) , and 2cr(3p) , 

Fig. 17. Ionization of carbon by electron impact. T(2p), T(2s), and T(2s + 2p) 

represent the present calculations and give cr(2p) , cr(2s) , and the sum 

of o-(2p) and cr(2s). No measurement is available on ionization of this 
1 1 

element 

Fig. 18. Ionization of atomic nitrogen by electron impact. Measurement of 

Smith et al is for total ionic current while that of Peterson is for 
% 

N current. For a description of the experimental results see 
Ref. 15. T(2p), T(2s), and T(2s + 2p) represent the present calcu- 
lations and givecr(2p), cr(2s), and the sum of <r(2p) andcr(2s'). P(2s + 

2p) is calculation of Peach (Ref. 4) for the sum of cr(2p) and'0'(2s) 
plus autoionization. 

Fig. 19. Ionization of atomic oxygen by electron Impact. Measurement of 

Fite et al and Bokesenberg is for O’*' production while that of Rothe 
et alls for total ion production. For a description see Ref, 15. 
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T(2p), T(2s), and T(2s + 2p) represent the present calculations and 
give cr{2p), cr(2s), and the sum of cr(2p) and cr(2s) respectively. S(2p) 
is calculation of Seaton (Ref. 37) for cr(2p). 


Fig. 20. Cross section per unit energy range of ejected electrons in ionization 
of carbon by electron impact. The energy of the primary is marked 
on each curve. 

Fig. 21. Cross section per unit energy range of ejected electrons in ionization 
of nitrogen by electron impact. The energy of the primary is marked 
on each curve. 

Fig. 22. Cross section per unit energy range of ejected electrons in ionization 
of oxygen by electron impact. The energy of the primary is marked 
on each curve. 

Fig. 23. Ionization of lithium ion by electron impact. The experimental cross 
section for production of Li‘‘'+ is due to the absolute measurement of 
Peart, Walton, and Bolder (Ref. 38), and Peart and Bolder (Ref. 39). 

R represents our calculation with relativistic correction. NR is our 
non-relativistic Born calculation. Presented also are calculations of 
Moores and Nussbaumer (Ref. 42) shown by MN curve in the figure. 
For other calculations see Ref. 38. 
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Appendix I 

Values of the Effective Charges Z e and the Ionization Potentials IP Employed^ 


Ltomic 

umber 

Element 

Shell 

Ze 

IP 

(eV) 

Atomic 

Number 

Element 

Shell 

Ze 

IP 

(eV) 

2 

He 

Is 

1.618 

24.580 

19 

K 

4s 

4.577 

4.339 








3p 

8.700 

24.63 

3 

Li 

2s 

1.549 

5.39 



3s 

10.571 

40.8 



Is 

2.617 

67.407 




, 







30 

Zn 

4s 

8,282 

9.391 

6 

C 

2p 

2.869 

11.264 



3d 

12.002 

17.4 

- 


2s 

3.784 

19.375 



3p 

17.366 

93.4 

7 

N 

2p 

3.456 

14.54 

36 

Kr 

4p 

11.785 

13.996 



2s 

4.524 

20.325 



4s 

14.730 

27.6 








3d 

19.06 

96.86 

8 

0 

2p 

4.035 

13.614 








2s 

5.260 

28.44 

37 

Rb 

5s 

6.659 

4.176 








4p 

13.257 

21.10 

10 

Ne 

2p 

5.180 

21,599 



4s 

16.009 

34.74 



2s 

6.726 

47.5 











54 

Xe 

5p 

15.612 

12.127 

11 

Na 

3s 

3.208 

5.138 



5s 

18.930 

23.4 


J 

2p 

6.262 

38.094 








2s 

7.702 

70.75 

55 

Cs 

6s 

8.564 

3.893 








5p 

17.289 

17.859 

12 

Mg 

3s 

4.150 

7.644 



5s 

20.387 

28.14 



2p 

7.299 

58.2 








2s 

8.691 

96.2 

3 

LI 

Is 

2.683 

75.282 

18 

Ar 

3p 

7.517 

15.755 








2s 

9.493 

29.24 







^The effective charge for atoms were taken from Ref. 12, and for Li'*' from 
Ref. 13. All outer shell ionization potentials were taken from M. A. Lange 
in the Handbook of Chemistry, 10th Ed. {McGraw Hill, Inc., New York, N.Y., 
1961) p. 111. The inner shell ionization potentials were taken from the 
available photo ionization data (Cf. R. D. Hudson and L. J. Kieffer, JILA 
Information Report No. 11, JILA, U. of Colorado, Boulder, Colorado, 1970). 

In other cases energy levels based on X-ray data were used (Cf. J. C. Slater, 
Phys. Rev. 1039 (1955), and J. A. Beardon and A. F. Burr, Rev. Mod. Phys. 

125 (1967)) 
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